Dephasing in two decoupled one-dimensional Bose-Einstein condensates and the 
subexponential decay of the interwell coherence 
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We provide a simple physical picture of the loss of coherence between two coherently split one- 
dimensional Bose-Einstein condensates. The source of the dephasing is identified with nonlinear 
corrections to the elementary excitation energies in either of the two independent condensates. We 
retrieve the result by Burkov, Lukin and Dernier [Phys. Rev. Lett. 98, 200404 (2007)] on the 
subexponential decay of the cocherence tx exp[— (t/to) 2 ^ 3 ] for the large time t, however, the scaling 
' of to differs. 

O ; 

■ PACS. 03.75.Gg - Entanglement and decoherence in Bose-Einstein condensates, 03.75.Kk - Dy- 
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Effectively one-dimensional systems of ultra cold atoms are a model systems to study the fundamental processes 
' of the coherent dynamics and (de) coherence in interacting many body systems. In addition in the limit of zero 
£NJ , temperature they are a primary example of the exactly integrable Lieb-Liniger model [H, 0] ■ 

Recently experimental progress on both, optical lattices [|| and atom chips 0, Q , allow to confine ultra cold atoms 
in strongly elongated traps with o> r 3> oj z (<-^>r, <-^>z being the frequencies of the radial and longitudinal confinement, 
• respectively). These traps are an ideal system for studying ID physics as long as both the temperature T and chemical 
potential fi are small compared to the energy scale given by the transverse confinement: /i < huj r , ksT < fku r . Optical 
lattices enable study of global properties of ensembles of Id systems down to very small atom numbers and into the 
strongly correlated regime, atom chips allow to study the properties and dynamics of single Id systems. 

Strong inhibition of thermalization, a signature of integrability, was observed with bosons deep in the ID regime 
@, and interference experiments on atom chips with pairs of weakly interacting Bose gases, easily fulfilling the above 
conditions for one dimensionality, allowed to study the dynamics of (de) coherence Q and the interplay between 
' thermal and quantum noise 0]. 

The special interest in the decoherence in such effectively one-dimensional ultracold atomic systems is rooted in 
I i the fact that they display decoherence and thermalization despite being at the first glance a prime example of the 
exactly integrable Lieb-Liniger model [HQ. In a recent paper [|| , we have shown that even if the temperature and 
chemical potential are well below the energy of the radial excitation (/z < huj r , ksT < huj r ) radial modes can be 
excited virtually. These virtually excited radial modes give rise to effective three-body velocity-changing collisions 
which lead to thermalization and the break down of integrability. The typical thermalization scale for typical Id atom 
chip experiments @, is in the order of 100 ms even when thermalization due to two body collisions is completely 

suppressed. However the dephasing dynamics observed in such coherently split ultracold Id atomic clouds is 
I ■ even one order of magnitude faster . 
| In addition, the experiment on the time evolution of interference between two coherently split one-dimensional (ID) 
atomic Bose-Einstein condensates (BECs) |7| revealed an surprising sub-exponential decay of the inter- well coherence 

O; (^(x,t)^(x,t))cxexp[-(^ ) Q ]. (1) 
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The measured decay exponent a sa 2/3 is in agreement with the theoretical calculations of Burkov et al. [12j which 
predicts a = 2/3. In their theoretical approach the inter- well coherence decay is treated in the terms of the heat flow 
between symmetric, ip + — [tfjr, + ^r]/^/2, and antisymmetric, ^_ = {ipr, + ^r]/V^, modes. The two individual (fully 
split) condensates are designated as right (R) and left (R), tPr^l(x, t) being the corresponding atomic field annihilation 
operator in the coordinate representation. 

In the present work we choose a different, more intuitive way of describing the system and consider dephasing of 
two independent integrable systems with inter-correlated initial conditions. As in Burkov et al. [HJ we consider the 
weak interaction limit, which, together with the finite size of the system, allows for a finite condensed fraction, in 
contrast to the Tonks limit or the case of an infinitely long quasicondensate. In what follows we use the system of 
units where Planck's and Boltzmann's constants are set to 1. 
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After the splitting, the system consists of two independent BECs, each being described by the Hamiltonian 
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where m is the atomic mass, gm is the ID coupling constant. The main contribution to the coherence reduction 
stems from phase fluctuation, so that 
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where the phase operators <f>± — ((f>Li(f>B,)/V2 and their conjugate density fluctuations operators n± — (ul ± ^i?)/\/2 
are defined via ipj = (n + fij) 1 ^ 2 exp(i(j>j), j = L, i?, n being the average ID BEC density. In what follows, we consider 
only classical (thermal-like) fluctuations and therefore omit operator notations, writing simply <pj etc. 
The local fluctuations can be expanded in plane waves, 
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C being the quantization length. The linearization of Eq. ^ yields Bogoliubov spectrum u k — 
yjk 2 j (2m)[k 2 j (2m) + 2mc 2 ], c being the speed of sound. In what follows we consider BEC at temperatures be- 
low the chemical potential, so that we assume phonon spectrum of excitations of the uniform ID Bose gas at rest, 
uj k = c\k\. 

The key idea of our treatment is to recall that the local density and velocity fields in BEC fluctuate, thus giving rise 
to the random non-linear corrections to the phonon frequency. Phonons propagating in the left and right condensates 
interact with different fluctuations, and this is the source of dephasing. 

The phonon energy depends on the BEC local density n = n + 6n via the speed of sound c = ygwn/m. Also the 
fluctuating local velocity V = m _1 90/ 'dx in the BEC contributes to the random energy correction as the advective 
term, so that the total correction to the average phonon energy c(n) \k\ is 
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Fluctuations of the density and local velocity differ for the right and left BECs, because the splitting process is never 
completely adiabatic (the degree of nonadiabaticity has been recently quantified by Polkovnikov and Gritsev [13]), and 
initial quantum (zero-point) fluctuations of (f>~ and n~ can be amplified to the values, comparable at low momenta 
to the initial thermal (classical) fluctuations of <fi + and n + . Still quantum fluctuations on their own can initiate the 
decay of the interwell coherence [l4| (on a different time scale). Since the density and local velocity fluctuations in 
the right and left BECs are different, the random energy shift is also different for the excitations propagating in these 
two BECs: 



Using Eq. Q, we obtain 



where 



fa(t) = exp( -ic\k\t-ij dt'SLu+) 0~ (0) cos J k (t) - i(p+(0) sin J k {t) 



J k (t)= / dt'8^{t'). 
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Assuming that fluctuations in fc-mode at t = are not correlated with the fluctuations of Su k at later times, we 
obtain from Eq. 
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In the short-time limit Eq. ([9]) reduces to 
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Moreover, since the system under consideration is closed and consists of two independent (L and R) integrable 
subsystems, the sum ((<j) k ) 2 ) + (((f>§ ) 2 = ({(f> k ) 2 ) + ((<P k ) 2 ) is time-independent, and we obtain 
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Since the correlation function for the random interwell frequency shift is 
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where 5lo£ characterizes the amplitude of the frequency shift fluctuations for the mode with given k, and /fe(r) is 
the correlation function in dimensionless form, equal to 1 at r — and rapidly approaching if its argument exceeds 
certain correlation time t£. Then Eq. (fTTj) yields 
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where the dephasing rate is 



Assuming that the condensate splitting is close to adiabatic and, hence, 

(fe(O)] 2 )«([0+(O)] 2 ) 
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and neglecting the change (approximately by a factor of 2) of ((fa) 2 
approximations are assumed in Ref. fl2l|). we obtain in the limit t > i 



in the course of the system evolution (similar 
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The remaining question is to determine Slu^ and r 
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at least roughly, we apply the following 



method. Instead of plain waves, we consider a wave packet centered at the momentum k and having the momentum 
uncertainty Ak ~ k. Such a width allows us to choose the shape of the wave packet close to the minimum- uncertainty 
wave packet, that allows us to localize it on a spatial scale of the order of 2ir jk. Then we can identify the random 
energy shift experienced by such a wave packet with that of the phonon with momentum k. 

In such a context, obviously, only the density and velocity fluctuations at wavelengths longer than 2ir/k contribute 
to 8lj^ . Fluctuations at shorter wavelengths are effectively averaged out and cause no influence to the dynamics of 
a wave packet of a spatial extension ~ 2n/k. Taking into account that fluctuations at different momenta are not 
correlated and relacing sum over discrete state by integration over continuous spectrum, we obtain 



8oj 



A 2 



2, 2 

z k 



k dk[ 
2tt 



-k 



((n- k >Y) , k 2 ((fa,f 



8n 2 



2m 



2/2 



(18) 



Here n^T is the Fourier component of the density fluctuations introduced analogously to Eq. . 

The remaining issue is to estimate the correlation time rf,. Say, the wavepacket propagates along x in the positive 
direction. Half of the surrounding fluctuations propagates in the opposite direction, bringing about a short correlation 
time scale ~ l/(cfc). However, half of the fluctuations co-propagate with the wave packet, in the first approximation 
at the same velocity. If there were no dephasing (with respect to each other) of fluctuations at different momenta, 
the corresponding correlation time would be infinite. However, longer-wavelength correlations dephase as well, at 
a rate similar to that of the wave packet under consideration. Averaging the dephasing rates over the ensemble 
of fluctuations restricted to \k'\ < \k\, we obtain the correlation time of the fluctuations affecting the dynamics of 
phonons with momentum k to be 

r c k « T- 1 . (19) 



The correlation time cannot be longer than given by Eq. (|19j) , because otherwise Eq. 
and the dephasing is slowed down significantly. Eqs. (fT51 [TT)|) result in 
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It is unlikely that deeper insight to the short-time dynamics of nonlierly interacting modes of a ID BEC can enhance 
the dephasing rate further compared to Eq. (f2"U|) . because short-time dephasing given by Eq. (fT3"|) is of the form 
exp(— const t 2 ) sa 1 — const t 2 , and fast perturbation can cause only slowdown of the evolution (Quantum Zeno effect) , 
but not speed up (anti-Zeno effect) [la ]. Eqs. (fTTI |2T)]) yield finally the following dephasing dynamics: 
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(the factor of 2 in the decrement is omitted, because of inexact nature of the estimations involved). 

In the asymptotically long time limit we may estimate fluctuation amplitudes from the ID Bogoliubov treatment 
for the phononic ensemble at temperature T as illustrated in [l6| (the same assumption is taken in Ref. [l2l]): 
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where [i — mc 2 is the zero-temperature approximation for the chemical potenial. Substituting Eq. (|22j) into Eq. (|18|) . 
we obtain 



(23) 



Finally, Eq. (|2Tj) takes the form 
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which is to a certain extent similar to Eq. (14) of Ref. (l2j ]: in both cases the de pha sing rate is proportional to |fc| 3//2 , 
but there is also an important difference: the relaxation rate predicted in Ref. [12| is by the factor of n/(mc) faster 
than our estimate. Eq. (f24|) has a generic form 
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where b\, 62 are certain constants and ZF{&) is a function that is finite at — > 00 and decreases faster than 0^2/3) if 
— > 0, to provide the convergence of the integral 
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that appears in the coherence factor Eq. The cut-off momentum k max depends on the temperature and chemical 
potential and is of the order of mc if T *~ /x. For t 3> T^ 1 , fi^ 1 we can substitute k max by 00. Changing the integration 
variable from momentum k to the dimensionless time 6 = b2\k\ 3 / 2 t we obtain 
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that provides the subexponential decay Eq. {1} with a — 2/3. Comparing Eqs. (|24|) and ([24| we find the scaling time 
to of the subexponential decay Eq. ([TJ) : 
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Note that the time limit t 3> 7 1 " 1 , /i -1 is not sufficient for Eq. (f2"Tl) to hold, since the latter is derived under assumption 
of significant thermalization of the antisymmetric mode that should happen at t ~ to. Our estimation of to is by 
the factor of 0.24 n/(mc) > 1 longer than that of Ref. [Hj] and the scaling with experimental parameters are also 
different. Recalling that fi — gujn with guj oc v±_ being the radial trapping frequency of the atomic waveguide 

and regarding it. the ID number density and experimentally obtained to as input parameters, we obtain T cx n/to, 
whereas Ref. 12] gives 

Let us now compare our results to the findings of Burkov, Lukin and Demler [121 ] in more detail: 
In the theoretical approach of Ref. [13], the inter- well coherence decay is treated in the terms of the heat flow 
between symmetric, ip+ = [ipL + ^r]/V^, and antisymmetric, ^_ = [i/>l + ^Jn]/y/2, modes. This point of view is 



counterintuitive for two completely split one-dimensional BECs which are two independent close to exactly inteqrable 
systems [l|, 0]. The thermalization times calculated, even when including the new virtual 3-body collisions Q, are 
much longer then the observed decoherence. The termalization time scale should be the one at which the heat flow 
arguments should not apply. 

The results of Ref. [l2| seem to imply some unphysical consequences: If we consider T ~ fi and estimate the 
dephasing rate for phonons with the energy close to the chemical potential Ref. [l2j gives a rate in the order 
of y/h/(mc)ij, 3> [i (since the experimentally accessible ID BECs are characterized by n/(mc) « 30), which is quite 
counterintuitive: the phonons with k ~ mc become over damped, and their dephasing rate exceeds any frequency scale 
available in the system with T ~ /x. Under the same condition our theory predicts dephasing rate ~ fi/ \Jn/ (mc) <C [i. 

The time scale to as suggested by our estimates implies the reconsideration of the temperature estimations for the 
experimental data of Ref. Q . Our model suggests that the actual final temperatures were higher by a factor of order 
2 than it was concluded from [T3|. This may easily happen, since the mechanisms of the heating of a BEC during 
the splitting (which is adaibatic only partially, taking into account its time scale ~ 10 ms) are not yet explored and 
understood. 

This difference between the prediction of the timescales and the scaling with experimental parameters between our 
model and the calculations by Burkov, Lukin and Demler [l2l ] demands both new more detailed experiments over a 
wider parameter range and an comprehensive numerical calculation of the splitting and coherence dynamics. New 
ways to measure temperature using the statistics of interference patterns 0] , and more refined RF potentials for atom 
manipulation [l7j will greatly extend the capability for experimental investigations. 
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